PDK-788 
TUHEP-02-01 
May 2002 



•l 



Search for Neutron- Antineutron Oscillations 
Using Multiprong Events in Soudan 2 



J. Chung 5 , W.W.M. Allison 3 , G.J. Alner 4 , D.S. Ayres 1 , W.L. Barrett 6 , 
P.M. Border 2 , J.H. Cobb 3 , H. Courant 2 , D.M. Demuth 2 
T.H. Fields 1 , H.R. Gallagher 5 , M.C. Goodman 1 , R. Gran 2 , 
T. Joffe-Minor 1 , T. Kafka 5 , S.M.S. Kasahara 2 , P.J. Litchfield 4 , 
W.A. Mann 5 , M.L. Marshak 2 , R.H. Milburn 5 , W.H. Miller 2 , L. Mualem 2 , 
A. Napier 5 , W.P Oliver 5 , G.F. Pearce 4 , E.A. Peterson 2 , D.A. Petyt 4 , 
K. Ruddick 2 , M. Sanchez 5 , J. Schneps 5 , A. Sousa 5 , 
B. Speakman 2 , J.L. Thron 1 , S.P Wakely 2 , N. West 3 

1 Argonne National Laboratory, Argonne, IL 60439 

2 University of Minnesota, Minneapolis, MN 55455 

3 Department of Physics, University of Oxford, Oxford 0X1 3RH, UK 
Rutherford Appleton Laboratory, Chilton, Didcot, Oxfordshire 0X11 OQX, UK 
5 Tufts University, Medford, MA 02155 
6 Western Washington University, Bellingham, WA 98225 

February 7, 2008 



Abstract 

We have searched for neutron-antineutron oscillations using the 5.56 fiducial kiloton-year ex- 
posure of the Soudan 2 iron tracking calorimeter. We require candidate nn occurrences to have 
> 4 prongs (tracks and showers) and to have kinematics compatible with TTN annihilation within 
a nucleus. We observe five candidate events, with an estimated background from atmospheric 
neutrino and cosmic ray induced events of 4.5 ± 1.2 events. Previous experiments with smaller ex- 
posures observed no candidates, with estimated background rates similar to this experiment. We 
set a lifetime lower limit at 90% CL for the nn oscillation time in iron: Ta(Fc) > 7.2 x 10 31 years. 
The corresponding lower limit for oscillation of free neutrons is T n n > 1-3 x 10 8 seconds. 



PACS numbers: 11.30.Fs, 12.20.Fv, 12.60.Jv, 14.20.Dh 



1 Introduction 



1.1 Neutron-antineutron oscillations 

An intriguing variation on the grand unification theme that nucleons are likely to be unstable 
is the proposal that neutrons can oscillate into antineutrons. Neutron-antineutron oscilla- 
tions were first predicted in 1970 by V. A. Kuzmin in a model intended as a realization 
of requirements given earlier by A.D. Sakharov for evolution of the universe to net baryon 
asymmetry ||. Subsequently nn oscillations emerged as a predicted reaction in certain 
grand unification theories ||. More recently it has been shown that nn oscillations can occur 
in a large class of supersymmetric SU{2) l x SU (2)# x SU (4) c models [4| . In such models the 
dominant baryon number violating process is a AB = —2, AL = nucleon transition (e.g. 
nn oscillations or p + n — > pions) rather than a AB = —1,AL = — 1 nucleon- antilepton 
transition (e.g. p — > e + ir or p — > VK + ). Neutron-antineutron oscillations have also been 
indicated as viable by recent GUT models which invoke the existence of extra spacetime 
dimensions ||. 

If indeed a neutron can evolve into an antineutron, the experimental signatures for the 
metamorphosis should be distinctive. The resulting antineutron will annihilate with a baryon 
of the surrounding environment, producing multiple mesons (B = 0) whose visible energy 
and net momentum are approximately those of two nucleon masses having nuclear Fermi 
motion. 

From the phenomenology of neutron-antineutron oscillations it can be shown || that the 
oscillation time T4 of a neutron bound within a nucleus of atomic mass A is related to the 
neutron oscillation time in vacuum r nS according to 



Here Tr, which has units of inverse time, is the suppression factor representing the effect of 
the nuclear environment which substantially prolongs the effective oscillation time. 

Detailed calculations of the suppression factor Tr for parent nuclei of experimental in- 
terest, including deuterium, oxygen, argon, and iron, have been reported in literature. The 
calculations utilize phenomenological frameworks provided by nuclear potential theory |7]. |J . 
and by S-matrix theory [H, [TU], |TT], |I2"|| . In the analysis of Dover, Gal, and Richard [0], it is 
proposed that neutron-antineutron oscillations will occur mostly in outer nuclear shells and 
near the nuclear surface. However reservations concerning this picture have been expressed 
and in a number of calculations the entire nuclear volume contributes to nn oscillations ||. 

1.2 Previous experimental searches 

Two types of experiments have been used to search for neutron-antineutron oscillations. 
In one approach, slow neutrons from a fission reactor are channeled through a magnetically 
shielded vacuum pipe towards a target region. An antineutron produced during the flight will 
annihilate in the target and the annihilation products are registered by detectors surrounding 
the target. Experiments of this type have been carried out at Pavia and at Grenoble 
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[L5[]. The Grenoble reactor experiment obtained r n7T > 0.86 x 10 8 s at 90% confidence 
level (CL). This is the most stringent oscillation time lower limit reported to date using free 
neutrons. 

The alternate approach, used in this experiment, is to continuously monitor neutrons 
bound in nuclei, usually as part of an ongoing nucleon decay search. Searches of this type 
have been reported by the underground experiments Homestake [TO, NUSEX [T7J], KOLAR 
18] , 1MB Kamiokande [20|, and Frejus [plfl . The searches by Kamiokande and Frejus 



obtained the most stringent nn oscillation time lower limits. 

In a search based upon a 1.11 kiloton-year (kty) exposure of the Kamiokande-I water 



Cherenkov detector p0[ , no candidate nN event was observed. An oscillation time lower 
limit of Ta > 4.3 xlO 31 years at 90% CL was set. Using the suppression factor Tr = 1 x 10 23 
s _1 calculated by Dover et al. for oxygen 0, Kamiokande obtained an oscillation time 
limit for free neutrons of r n n > 1.2 x 10 8 seconds at 90% CL. The Frejus collaboration, in a 
search using a 1.56 fiducial kty exposure of the experiment's planar iron tracking calorimeter, 
also reported zero nn oscillation candidates. The oscillation time lower limit for Ta in iron 
thereby obtained was 6.5 x 10 years at 90% CL. Using T R = 1.4 x 10 23 s" 1 as calculated 
by Dover et al. for iron, Frejus determined the free neutron limit to be also r n7 j > 1.2 x 10 8 
seconds at 90% CL 

In the Kamiokande analysis an enhanced probability for nn oscillations to occur in the 
nuclear periphery as postulated by Dover et al. was assumed. In the Monte Carlo simulations 
of the experiment, the effect of this assumption is to reduce distortion of the final state meson 
spectrum arising from intranuclear absorption and inelastic scattering processes. As a result 
experimental detection efficiencies are enhanced relative to expectations for the case where 
oscillations may occur throughout the entire volume of parent nuclei. The search reported 
here follows the more conservative approach adopted previously by Frejus. For our primary 
nN simulation on which our detection efficiency is based, we assume nn oscillation to occur 
throughout the nuclear volume. However the effects of peripheral predominance are also 
described. 



2 Detecting nn Oscillations in Soudan 2 

2.1 Detector and data exposure 

Soudan 2 is a 963 metric ton (770 tons fiducial) iron tracking calorimeter of honeycomb 



lattice geometry which operates as a slow-drift time projection chamber [p^j . The tracking 
elements are one-meter-long drift tubes filled with an Argon-C02 gas mixture. Electrons 
liberated by throughgoing charged particles drift to the tube ends under the action of a 
voltage gradient applied along the tubes. The drift charge is registered by vertical anode 
wires, while horizontal cathode pad strips register the image charges. The third coordinate 
is obtained from the drift-time. The amount of charge measures the deposited ionization. 
The drift tubes are laid onto corrugated steel sheets, and the sheets are stacked to form 
1x1x2.5 m, 4.3 ton modules from which the calorimeter is assembled in building-block 
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fashion. Surrounding the tracking calorimeter on all sides but mounted on the cavern walls 
and well separated from calorimeter surfaces is a 1700 m 2 active shield array of two or three 



layers of proportional tubes | Z3|| . The shield facilitates identification of events which are not 
contained within the calorimeter. In particular, it provides tagging of background events 
initiated by cosmic-ray-induced neutrons. 

The detector is located at a depth of 2070 meters-water-equivalent on the 27th level of the 
Soudan Underground Mine State Park in northern Minnesota. The modular design enabled 
data taking to commence in April 1989 when the detector was one-quarter of its final size; 
routine operation with the fully-deployed detector got underway in November 1993. The 
fiducial (total) exposure analyzed here, obtained from data-taking through December 2000, 
is 5.56 (6.96) kiloton-years. 

Calibration of calorimeter module response was carried out at the Rutherford ISIS spal- 
lation neutron facility using test beams of positive and negative pions, electrons, muons, and 
protons |24] . Spatial resolutions for track reconstruction and for vertex placement in anode, 
cathode, and drift time coordinates are of the same scale as the drift tube radii, ~ 0.7 cm. 
In Soudan 2, ionizing particles having non-relativistic as well as relativistic momenta are 
imaged with dE/dx sampling in a fine-grained honeycomb lattice geometry. Protons can be 
distinguished from pions and muons via ionization and ranging, energetic muons discrimi- 
nated from pions via absence of secondary scatters, and prompt e^ 1 showers distinguished 
from photon showers on basis of proximity to primary vertices. These event imaging ca- 
pabilities offer advantages, in comparison to water Cherenkov detectors and to planar iron 
calorimeters, for analysis of the complicated multiprong topologies that would arise with nN 
annihilations initiated by nn oscillations. 

2.2 Simulation of nN annihilation arising from nn oscillations 

We have developed realistic simulations of nn oscillations yielding nN annihilations as they 
would occur in the iron nuclei which comprise the bulk of the calorimeter mass. Generation 
of nN events is carried out as follows: 

Momenta of the initial state antineutron and nucleon are assigned according to a dis- 
tribution based upon a Fermi-gas model parameterization of quasi-elastic electron nucleus 
scattering ||25|| . Final state particle four-momenta were constructed in accordance with N- 
body phase space [E3]. Assignment of nN reactions to generated events is weighted according 
to cross section data for pp annihilation at rest E7I, EH] . For the purpose of channel selection, 



nn annihilations are assigned the same cross sections as observed for pp; np is assumed to 
have the same total rate as pp, and cross sections were inferred from pp as allowed by charge 
conservation. We restricted our reaction compilation to cross sections exceeding 2% of the 
total NN cross section. Consequently production of p, to, charged and neutral pions and 
kaons is represented, however final states with rj and rj' mesons were neglected. 

Provision was made to include intranuclear rescattering (INS) of final state pions in the 
simulations. Our treatment follows the approach utilized previously in simulations of atmo- 
spheric neutrinos and of nucleon decay [Q. Pions, either directly produced in annihilation 
or in p decays, are propagated through a model nucleus. The nuclear radius and density are 
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parameters of the model; the radius is scaled according to A 1 / 3 . Scattering in the nuclear 
medium is characterized using a momentum-dependent pion interaction length. The nuclear 
parameters were set by requiring the model to reproduce single and multiple pion produc- 
tion rates observed in bubble chamber z/^-deuteron (A=2) and z/^-neon (A=20) interactions 
p9| . For pions from vertices placed at random within an iron nucleus, our model predicts 
~ 45% to emerge either unscattered or to have undergone only small angle elastic scattering. 
About 30% are predicted to undergo total absorption, while the remainder undergo charge 
exchange (9%) or emerge having undergone inelastic scattering (16%). Our treatment of 
nuclear rescattering does not include proton secondaries. While it is expected that protons 
would occasionally be ejected from parent nuclei as the result of rescattering, their momenta 
would almost always fall below the effective threshold (approx. 450 MeV/c) for creation of 
distinct tracks in the detector. 

For simulations with intranuclear rescattering included, the average multiplicity per event 
for mesons emerging from an iron nucleus is 3.8 with rms deviation 1.2. The average track 
plus shower multiplicity emitted from an iron nucleus ("prong multiplicity") per event is 5.0 
with an rms deviation of 2.1 prongs. (Note these multiplicities are for "perfect detection".) 

The nN event generation produces four-momenta of all particles which exit the parent 
nuclei. This information is fed, event-by-event, into the Soudan 2 Monte Carlo (MC) which 
provides a realistic detector response. The simulation includes detector background "noise" 
arising from natural radioactivity and from the electronics. Event records are generated with 
format identical to that of data events, allowing nn oscillation events to be processed using 
the same codes and procedures as for data and for events of the atmospheric neutrino Monte 
Carlo. 

2.3 Properties of processed nN samples 

Three separate nN simulation samples were generated. For each sample, events were gener- 
ated at random locations throughout the tracking calorimeter as it evolved during the ex- 
periment's eleven years of data-taking. In our "primary" simulation, nN annihilations were 
started at random points throughout their parent nuclei and pion intranuclear rescattering 
was implemented. In a second simulation which also included intranuclear rescattering, nN 
annihilations were restricted to the nuclear periphery (0.75-R < r < R, where R is the nu- 
clear radius). A third simulation was carried out in which no pion intranuclear rescattering 
was included. 

Each of the three oscillation samples was processed using a sequence of selections and 
procedures very similar to that routinely used in reduction of data events in Soudan 2. 
Each sample was subjected (in software) to the hardware trigger. Events which passed the 
hardware trigger requirements were subjected to two different software "Filter" codes. The 
Filters impose event containment criteria, e.g. that no track from the event approaches closer 
than 20 cm to the detector outer surfaces; they also mitigate against backgrounds arising 
from cosmic ray muons, natural radioactivity, and detector noise. Events which survived the 
Filters were then subjected to two separate scans by physicists. Scanning was carried out 
using interactive color graphics workstations. The scan rules provided refinements to the 
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Filter selections and introduced requirements on imaging quality, e.g. an event was rejected 
if (?) its primary vertex occurred in material interior to the detector but not instrumented, 
or if (ii) its proximity to inter- module gaps compromised the reconstruction of the event. 

Events which survived successive application of the hardware trigger, software filters, and 
physicists' scans are tallied in the upper four rows of Table p] for each of the three simulations. 
These entries are input to the calculation of the detection efficiency for neutron-antineutron 
oscillations in Soudan 2 as will be described in Sect. 



Event Processing and 
Selection Stages 


nN annihilation 
without INS 


nN annihilation in 
nuclear periphery 
with INS 


nN annihilation in 
nuclear volume 
with INS 


Initial sample 


491 


491 


491 


Hardware trigger 


490 


469 


451 


Containment and quality niters 


288 


301 


286 


Two physicist scans 


214 


229 


205 


Multiplicity > 4; no proton events 


172 


148 


135 


Exclude events with "muons" 


146 


130 


123 


Kinematic selection on E V i S , P/E 


137 


102 


86 



Table 1: Survival through successive processing and selection stages for events of three different simulations 
of nn oscillations yielding nN annihilations in Soudan 2. 

Preliminary to kinematic reconstruction, the topology of each event was characterized in 
terms of track and shower prongs. A track prong results when an ionizing, non-showering 
charged particle (e.g. a pion, muon, or proton) traverses drift tubes of the calorimeter's 
honeycomb lattice leaving a continuous trail of tube "hits" . A shower prong on the other hand 
is created by photon conversion or by a primary electron or positron. The electromagnetic 
shower consists of many distinct particles and exhibits a cone-like pattern of hits, interspersed 
with gaps, which is aligned with the direction of the initiating photon or electron. The mean 
conversion length for photons in Soudan 2 is approximately 15 cm. As a consequence, 
photon-induced showers generally appear in the vicinity of but moderately displaced from 
event primary vertices. 

Images of two nn oscillation events simulated with full detector response are displayed in 
Fig. 0, where the anode (X) versus drift time (Z) view of each event is shown. For event 
scanning, three views are always used, including cathode-time and anode-cathode as well as 
anode-time projections. (The anode-cathode view is automatically assembled for all events, 
MC as well as data, using demultiplexing and hit-matching algorithms.) As suggested by 
Fig. |I|, nN events appear to be energetic yet isotropic to a degree uncommon for data events 
and for events of the atmospheric neutrino Monte Carlo. 

The occurrence of relatively high multiplicity per event of track plus shower prongs is a 
signature feature of nN annihilation. From simulation of nn oscillation events in the absence 
of intranuclear rescattering, we find the prong multiplicity distribution - after triggering, 
filtering, scanning, and reconstruction - to have a mean of 6.3 prongs per event with an rms 
width of 2.8 prongs. Here the prong multiplicity is relatively high (compared to multiplicities 
as generated, see Sect. 2.2), due to depletion of low multiplicities by the processing. With 
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Figure 1 : Anode versus drift time views (magnified) of two simulated annihilation events within iron nuclei 
following nn oscillation. The multiprong hadronic final states of these events are isotropic to a degree 
unusual for neutrino interactions. 



intranuclear rescattering included, the prong multiplicity distributions of processed samples 
are shifted lower. For both simulations with INS included, the multiplicity distributions have 
mean values of 4.8 prongs with rms widths of 3.2 prongs. Thus a sizable contribution from 
multiplicities exceeding three-prongs survives in the latter simulations. Events of this kind 
are relatively uncommon among contained atmospheric neutrino events. 

In addition to event topology, the reconstructed kinematic quantities E vis , the visible 
energy, and P ne t, the net momentum, provide discrimination between nN and background 
events. These quantities and their correlation have been used previously in our searches 
for nucleon decay [pPj] . E vis is calculated as the sum of the relativistic energies for each of 
the final state tracks (using the pion mass) and showers. There occur a few short, heavily 



ionizing tracks which satisfy our identification criteria for protons ||31|| . For these, only the 
kinetic energy is added into the calculation of E vis . Distributions of E vis versus P net for nN 
events in the absence of and including intranuclear rescattering, are plotted in Fig. |2|a and 
Fig. |2|b respectively. Reconstructed events, even in the absence of INS as in Fig. |2]a, exhibit 
large energy losses arising from low energy prongs and secondary hadronic scatters which 
are unresolved. Consequently the events cluster well below E vis of 1.88 GeV. Comparison of 
Fig. ^|a with Fig. Qb shows that a further large degradation arises with INS. The INS effects 
substantially increase the kinematic overlap of nN events with atmospheric neutrino events. 



3 Data and Backgrounds for the nN Search 
3.1 Neutrino, rock, and atmospheric v MC events 

Since nN events generally have high multiplicities we restrict our analysis to events having 
a "multiprong" topology, that is, having two or more produced particles emerging from 
primary vertices. Quasi-elastic neutrino interactions which produce a charged lepton plus 
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Figure 2: Distributions of visible energy versus magnitude of vector net momentum for reconstructed nN 
annihilation events of two independent full detector simulations. The simulation plotted in |^a (||b) omits 
(includes) intranuclear rescattering of final state pions. The shift to lower E V i S in simulation of b) versus a) 
arises from pion absorption and inelastic scattering within iron nuclei. 



a recoil proton are readily distinguished and removed. As was done for the simulated nn 
events, we require all events under consideration to be fully contained in an interior volume 
which is everywhere 20 cm from calorimeter outer surfaces. 

Three multiprong event samples, which are in addition to the simulated nn samples dis- 
cussed above, have been isolated for this search: 

First and foremost are data events for which the cavern-liner active shield array was qui- 
escent during the allowed time window. These events comprise our shield-quiet data sample. 
They originate mostly with atmospheric neutrino interactions but would also include nn 
oscillation events. There are 188 multiprong events in the shield-quiet data of this analysis. 

There is a second category of data events which comprises a background to neutrino 
multiprongs as well as to nn oscillations. Events of this category are usually produced by 
energetic neutrons released in inelastic cosmic ray muon interactions with the cavern rock 
surrounding the tracking calorimeter. Most of these rock events are accompanied by charged 
particles which give coincident hits in the active shield array; they constitute our shield-tagged 
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rock event sample. However, a few rock events are not accompanied by shield hits, either 
because of shield inefficiency (detection efficiency is 94%) or because no charged particles 
enter the cavern along with a rock neutron. The latter shield- quiet rock events may end 
up as background in shield-quiet data. Fortunately the shield-tagged rock events provide a 
control sample from which the amount of residual rock background in shield-quiet data may 
be inferred. 

The predominant background for nn oscillations arises from the atmospheric neutrino 
reactions. Their contribution is evaluated using the Soudan 2 atmospheric neutrino Monte 
Carlo sample. Our neutrino MC simulation utilizes the atmospheric flux calculation of 



the Bartol group for the Soudan site [32]. Details of the neutrino event generation and 



comparisons with low energy uN data are given elsewhere [33 . MC atmospheric neutrino 



events were inserted into the experiment's data stream during data-taking. The MC events 
were introduced at a rate 6.06 times higher than the data rate expected (for null neutrino 
oscillations) from the atmospheric neutrino flux. The events were processed identically to 
real data events, with their identity revealed only at the final analysis stage. The atmospheric 
v MC multiprong sample finally extracted contains 1267 events. 

Properties of our multiprong data, rock, and neutrino MC samples germane to an nn 
oscillation search are reported below. Further details can be found in previous publications 
VA 1301, BH. 



3.2 Multiprong data compared to atmospheric v MC 

The prong multiplicity distribution of the shield-quiet data sample is summarized using tallies 
in the grid displayed in Table 0. In the Table, the topology of each event is represented by 
the number of tracks, n tra ck, and number of showers, n s hower, which comprise the visible final 
state. The number of data events which have a particular (n tr ack, n s hower) combination are 
given by the upper entries at the grid location having the appropriate integer coordinates. 
The multiprong data topologies having highest rates are seen to be combinations of lowest 
multiplicity, namely (n track ,n shower ) = (2,0) and (1,1). For comparison, the corresponding 
prong distribution for the atmospheric neutrino MC event sample is given by the lower 
entries in each of the track-shower grid locations of Table 0. The MC sample (for null 
neutrino oscillations) has been normalized to the same fiducial exposure as for the data 
(1267/6.06 = 209.1 events). The two distributions are seen to match well, e.g. the x 2 /bin 
averages to below 1.0, with the exception of the (1,2) topology. Since most shield-quiet 
data multiprongs are initiated by atmospheric neutrinos, the agreement is evidence that the 
neutrino MC provides a good general representation of our data. 

Kinematics for the shield-quiet multiprong data is shown by the E vis versus vector net 
momentum diplot of Fig. |]a. The events populate a broad, correlated region extending from 
threshold to 1.5 GeV (1.5 GeV/c) in visible energy (net momentum). Portions of this data 
distribution overlap the distribution predicted for nn events shown in Fig. |2|b. 

A direct comparison with kinematics for the atmospheric neutrino MC sample is provided 
by Fig. Pd. Note that the latter sample has an exposure-equivalent of 33.7 fiducial kton- 
years. The distribution of the MC sample, normalized to the exposure, is very similar to 
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Table 2: Event topology distributions: Event counts are tabulated according to (ntrack, ^shower) combination 
per event, for multiprong events of the data sample (upper entries) and of the atmospheric neutrino MC 
sample (lower entries) . Event tallies of the latter sample represent expectations for an atmospheric flux with 
null oscillations normalized (factor of 6.06) to the data exposure. 



that of the data. 

In Fig. ^, events which have prompt protons are denoted by solid circles. As previously 
noted, such events are highly improbable as nn oscillations. However our atmospheric v MC 
indicates that detectable recoil protons are to be expected in 24% of neutrino multiprongs, 
and this expectation is born out by the frequency of recoil proton events observed in the 
data of Fig. |3|a. In Fig. ^|a there is a tendency for data events having recoil protons to 
occur in a region parallel to and below the more populated kinematic band. This trend is 
reproduced in the MC distribution of Fig. |3|b. 



3.3 Background from rock events 

Rock data events which are flagged by hits in the surrounding shield array, have a distribu- 
tion of vertex depth into the detector which reflects the mean hadronic interaction length of 
neutrons in the Soudan detector medium, approximately 80 cm. A rock background compo- 
nent in our shield-quiet multiprong data will have a corresponding vertex depth distribution 
in the detector. Neutrino and nn oscillation events on the other hand, will distribute uni- 
formly throughout the detector volume. Thus by fitting the distribution of vertex depths 
observed in our multiprong data to the sum of two component distributions representing ver- 
tex depths of shield-tagged rock and of neutrino MC events, the amount of shield-quiet rock 



contamination can be determined IBTfl. On the basis of a maximum likelihood fit, we estimate 
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Figure 3: E V i S versus P ne t distributions for shield-quiet data events and for events of the atmospheric neutrino 
MC sample. Events without (with) proton tracks are depicted using open (solid) circles. The v MC sample 
(||a), which corresponds to an exposure (for null v oscillatiions) of 6.06 times the data, reproduces the general 
kinematic trends of multiprong data (^b). 



that a fraction 0.06to; 06 °f shield-quiet data multiprongs are produced by rock neutrons. 
4 Selection of nn Oscillation Events 

Four selection criteria for neutron oscillation candidates (Table ^) have been defined and 
applied to all event samples. The selections are designed to maximize nn detection in the 
shield-quiet multiprong data sample while minimizing backgrounds, the principal one arising 
from inelastic interactions of atmospheric neutrinos. 



Cut 


Definition 


1. 


N prong ^ 4 


2. 


No prompt "proton" tracks 


3. 


No prompt "muon" tracks longer than 150 cm 


4. 


Evt kinematics: P/E < 0.6, 700 < E ms < 1800 MeV 



Table 3: Selection cuts for nn oscillation events in Soudan 2. 



11 



1. Prong multiplicity: We require the sum of track and shower prongs directly asso- 
ciated with the primary vertices to be greater than or equal to four. This value is 
below but near to the mean multiplicities for processed nn samples as discussed 
in Sect. \2.'3\ Lowering our selection to include three-prongs would increase the 
efficiency for nN detection by only 3%, however backgrounds arising from atmo- 
spheric neutrinos and rock events would increase by 21% and 32% respectively. 

2. Primary proton tracks: Events having proton tracks emerging from their primary 
vertices are removed from consideration as nn oscillation candidates. The rate 
of prompt protons is negligible for nN annihilation events, however ~ 24% of 
atmospheric neutrino events and ~ 33% of rock events have visible recoil protons. 

3. Primary muon tracks: Pions produced in nN annihilations usually scatter, be- 
come absorbed, or range-to-stopping over distances comparable to the calorime- 
ter's hadronic interaction length of 80 cm. In order to mitigate against charged 
current background, we regard any non-scattering track which has pion/muon ion- 
ization and range-to-stopping > 150 cm (240 g/cm 2 ) to be a muon track. Any 
event having such a track is rejected. This selection eliminates 36% of the atmo- 
spheric neutrino sample while cutting less than 10% of nn events of the primary 
simulation. 

4. Event kinematics: As our final criterion we require the kinematics of nn candidates 
to be compatible with those of nN annihilation. The event net momentum fraction 
P ne t/E vis (hereafter P/E) has previously been shown to be a useful variable for 



separating nn oscillation events from atmospheric neutrino reactions plj , and its 
utility is confirmed in our analysis. Figures |]a and |]b show Pj E versus E V i S for nN 
events after topology cuts (1 through 3) from simulations in which oscillations occur 
throughout the parent nuclear volumes. For the case of no intranuclear rescattering 
(Fig. ||]a), approximately 94% of nn oscillation events occur in a region having Pj E 
less than 0.6, with E vis ranging from one to two nucleon masses. With intranuclear 
rescattering included, Fig. |]b, the nn event population shifts to lower values of 
E vis and to higher values of Pj E. Nevertheless a degree of clustering remains which 
allows a useful search region in the (E vis , P/E) plane to be defined. Our selection 
for nn kinematics is depicted by the rectangular region (dotted-line boundary) 
shown in the Figs |]-|6|. Candidate nn events are required to have P/E < 0.6 and 
700 < E V i S < 1800 MeV. These cuts are designed to minimize the atmospheric 
neutrino background shown in Fig. [5] while maximizing the acceptance for nN 
events. 



5 nn Candidates and Backgrounds 

The reductions in each nN annihilation sample upon successive nn candidate selections are 
summarized in the bottom three rows of Table [l|. It can be seen that intranuclear rescattering 
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Figure 4: Momentum fraction versus visible energy for events from nn simulation after topology cuts, with 
intranuclear rescattering omitted (Qa) versus included (jijb) . The kinematic search region of this experiment 
enclosed by the dotted-line boundary, contains 94% of events which pass topology selections. In the more 
realistic simulation of ||b however, INS degrades E V i S and increases momentum anisotropy on average, leaving 
70% of events in the search region. 



significantly lowers the survival rates. For our primary simulation, with INS operative and 
with annihilations occurring throughout the volume of the parent nucleus, 86 events from 
an initial 491 event sample survive all selections, giving a detection efficiency of 0.18 ± 0.02. 
Detection efficiencies for the comparison simulations are higher, (28 ±3)% for the simulation 
without INS, and (21 ±2)% for the simulation which includes INS but restricts annihilations 
to the nuclear periphery. 

Shield-quiet data events which satisfy nn selections 1-3 of Table |3] distribute in the P/ E 
versus E vis plane as shown in Fig. |6|. Of these sixteen events, five occur within the kine- 
matically allowed region and are candidate nn oscillation events of this search. Three of the 
candidates have four-prong topologies, while the remaining two are five-prong events. 

Two views of one of the candidates consisting of four pion-or-muon tracks emerging cleanly 
from a reaction vertex, are shown in Fig. [7]. This event is also consistent with multipion 
production by atmospheric neutrinos, e.g. © + n — > /i ± 7r T 7r + 7r _ (A r ) MM. 

Fig. [5] shows that, among atmospheric v MC events which satisfy topology selections 
1-3 of Table ||], there are 28 events (11 CC, 7 v e CC, and 10 NC) which occur in the 
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Figure 5: Kinematic distribution of atmospheric neutrino MC events which survive topology selections for 
nn candidates. Events which occur within the bordered region represent the background expected from an 
exposure (assuming null neutrino oscillations) which is 6.06 times the data exposure. 
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Figure 6: Distribution of data multiprong events after topology selections, in the P/E V i S versus E V i s plane. 
Five events are observed to have kinematics compatible with nn oscillations in the Soudan 2 detector. 



selected kinematic region. That is, 28 neutrino-induced background events would occur in 
a 33.7 fiducial kiloton-year exposure to an atmospheric neutrino flux having null neutrino 



oscillations. However strong evidence for atmospheric neutrino oscillations now exists [[35 
36fl , [31] , |37jj . We allow for neutrino oscillations in our background estimate by weighting 



our MC CC prediction by the atmospheric neutrino flavor ratio-of-ratios: R v = [(z/ M + 
Vfj) I {y e + V e )}E> a ta/[(vn + F M ) / [y e + V e )]Mc- We have used the value of R v measured in this 
experiment, based upon a 5.14 fiducial kton-year exposure: R v = 0.68±0.11±0.06 J37]. This 
correction reduces the neutrino background sample to 24.5 events. After scaling to the data 
exposure and accounting for statistical error plus uncertainties arising from our atmospheric 
neutrino MC including the neutrino oscillation correction (see Sect. |7.2| ), we obtain 4.0 ±1.0 
events as the estimated neutrino-induced background of our nn oscillation search. 

As discussed in Sect. |3.1| , neutron-induced rock events which are background for nn os- 
cillations are those which elude tagging by the active shield. As described in Sect. |T3 
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Figure 7: Anode-time and cathode-time views of a four-track nn oscillation candidate. 



we estimate that a fraction 0.06to;o<3 °f data multiprongs are rock events. Thus of the 188 
shield-quiet data multiprongs, we estimate ll.St^J to be rock events. However, among the 
375 shield-tagged rock multiprong events, only fifteen pass our four nn selection criteria. 
Thus, assuming that the selection efficiency of shield-quiet rock events is the same as that 
of shield-tagged rock events, the background of rock events which pass all nn selections is 
estimated to be 0.5±oJ events. 

In summary, we estimate the number of background events to be 4.0 ±1.0 (atm. v) plus 
0.5±g;l (rock) = 4.5 ± 1.2 events. 



6 Neutron Oscillation Time Lower Limits 

Having observed five candidate nn oscillation events and having calculated background neu- 
trino and cosmic ray processes to contribute 4.5 events, we set a lower limit for the nn 
oscillation time T4 in iron nuclei. We follow previous experiments in using our best estimate 
of the background in determining the limit. Additional uncertainties in the limit due to the 
uncertainties in the background calculation are discussed in Sect. |7|. We use the approach 
to confidence level construction as formulated by Feldman and Cousins PSL 1311. For five 



candidate events with 4.5 background events the signal limit at 90% CL is ngo = 5.5 events. 
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nn Simulation 


T4 lower limit 


r n n lower limit 


nn without INS 


11. x 10 31 y 


1.6 x 10 s s 


nn in Nuclear Periphery, with INS 


8.4 x 10 31 y 


1.4 x 10 s s 


nn throughout Nucleus, with INS 


7.2 x 10 3i y 


1.3 x 10 s s 



Table 4: Neutron-antineutron oscillation time lower limits (90% CL) arising from three different treatments 
of pion intranuclear scattering in nN annihilations. Simulations of the second and third rows are regarded 
to be realistic; the one with more conservative nuclear modelling (bottom row) is the basis for the limits for 
this search. 



The oscillation time lower limit is then calculated using 

T A > (2) 

Here, N n = 3.15 x 10 32 neutrons in a kiloton of the Soudan 2 detector, Tf = 6.96 kiloton 
years is the full detector exposure, and e = 0.18 is the efficiency for detection of nn oscillations 
in the full detector mass, calculated using our primary nn simulation (see Sect. ^). At 90% 
CL we obtain T A > 7.2 x 10 31 years. 

Equation ([I]) shows that the free nn oscillation time r n ^ can be inferred from oscillation 
time T4 for neutrons of nucleus A. For the suppression factor of iron we use the value 
Tr = 1.4 x 10 23 s _1 obtained by Dover, Gal, and Richard upon averaging over calculations 
using different nuclear optical potentials [0]. Similar or moderately higher Tr values have 
been reported from other calculations || [LI]. At 90% CL we obtain r n ^ > 1.3 x 10 8 seconds. 

Oscillation time lower limits based upon detection efficiencies of each of our three nn 
simulations are summarized in Table f|. The range of values which results from the variation 
in e illustrates the important role of the intranuclear rescattering treatment. The alterna- 
tive simulations serve as estimators of sensitivity to systematic uncertainties arising from 
modelling of the nuclear environment. 



7 Statistical and Systematic Uncertainties 

Our lower limits for Ta and r n ^ contain uncertainties arising from finite statistics of data 
and simulation samples and from systematic errors inherent in our analysis techniques. To 
quantify these uncertainties, we have evaluated the error contributions which enter via each 
of the factors of Eq. @ . 



7.1 Exposure and nn detection efficiency 



The N n value [^j an d the detector exposure are known accurately, with a consequent frac- 
tional error for the factor N n x Tf of order 2%. 

The nn detection efficiency e has statistical error due to the finite sample size of our 
primary simulation at generation and especially after cuts. Furthermore, systematic un- 
certainties may arise due to inaccuracies in the simulation. Based upon trial variations of 
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relative rates among dominant nN channels (cross sections known to < 20%), we estimate 
channel rate uncertainties to contribute ~ 5% uncertainty to e. A systematic error in e may 
arise from our placement of annihilation sites within parent nuclei. Based upon our alter- 
native simulations, we assign a 17% fractional uncertainty to modelling details. There are 
uncertainties inherent with our treatment of pion intranuclear rescattering. Compared to the 
unphysical case wherein INS is neglected, our treatment gives rise to an efficiency decrease 
Ae = 0.10. Uncertainty in the amount of this decrease arises from finite statistics of bubble 
chamber neutrino samples |HJ and with extrapolation from A = 20 to 56 of the component 
INS processes, namely pion absorption, charge exchange, inelastic and elastic scattering. We 
estimate the latter uncertainties to total 30% |HJ . We then infer the fractional error contri- 



bution to e arising from our INS treatment to be 17%. Combining the above uncertainties 
in quadrature, we estimate the total fractional uncertainty 5e/e to be 27%. 



7.2 Background estimation 

The value of ngo relies upon our estimation of rates for background events which satisfy the 
nn oscillation search criteria. 

Our estimate of atmospheric neutrino backgrounds is susceptible to errors from three 
sources: (i) Uncertainty arises from event statistics (after cuts) of the atmospheric neutrino 
simulation (19%); (ii) The normalization of the v MC to the experimental exposure, which 
predicts 183 neutrino events will be observed, is uncertain. From our data we observe (data- 
rock) = 177 d= 20 events, from which we infer an error (11%). (Hi) The atmospheric v MC 
may not fully represent aspects of neutrino data which feature in the selections of Sect. f|. 
However our data imply limits on the extent of MC mis-representation. For data events, 
and with respect to the cuts applied in Sect. f|, we observe 34 ± 5% to have n prongs > 
4, 78 ± 9% to be devoid of protons, 67 ± 8% to be devoid of "muons", and 43 ± 6% 
to have kinematics in the vicinity (500 < E vis < 1800 MeV, P/E < 0.8) of our search 
region. The corresponding fractions for the v MC sample are 35%, 76%, 64%, and 39% 
respectively. The agreement is good, and we take the quadrature sum of the fractional 
differences | Data — MC \ /Data in the four ratios as our estimator of uncertainty arising 
from MC representation of atmospheric neutrino physics (11%). The total fractional error 
for our atmospheric neutrino background is then 25%, which corresponds to an error of ± 
1.0 events on our estimate of 4.0 neutrino-induced background events. 

Concerning our estimate of background from rock events, the uncertainty arising from 
the determination of rock-event contribution to the vertex depth distribution of shield-quiet 
multiprongs outweighs any other systematic uncertainty. 

The total absolute error assigned to our 4.5 background events from atmospheric neutrinos 
and rock amounts to ± 1.2 events. Variation of our background estimate by ±1 a yields 
5ngo/n$Q = 21%. The total uncertainty in our oscillation time lower limit value for Ta 
thereby implied by Eq. (fj) is 5Ta /Ta = 34%. 

Our limit for the oscillation time of free neutrons r nS depends upon the nuclear suppres- 
sion factor Tr as implied by Eq. (0). For Tr of iron, we used the value of Dover et al. J7|, 
Tr = 1.4 x 10 23 s -1 . However from the range of Tr (A=56) values indicated by the calcula- 
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Experiments 


Source of 
neutrons 


Exposure 
fiducial (total) 
kton-years 


Cand. 
Events 


Est. 
Bkgrd 


Tr 
10 23 s" 1 


T A 
10 31 y 


7~nn 

10 8 s 


Grenoble('90) [15| 


reactor beam 














0.86 


Kamiokande('86) @ 


16Q 


1.11 





1.2 


1 


4.3 


1.2 


Frejus('90) |l| 


56 Fg 


1.56 (2.56) 





2.5,2.1 


1.4 


6.5 


1.2 


Soudan 2 [this study] 


5B Fe 


5.56 (6.96) 


5 


4.5 


1.4 


7.2 


1.3 



Table 5: Experimental lower limits at 90% CL for nn oscillation times Ta and T n n of bound and free 
neutrons respectively. For results of this search (bottom row), the background estimate is corrected for 
Vfj, — > v T oscillations and limits are calculated using the Feldman- Cousins method. 



tion of Alberico et al. [[Tl|], a theoretical uncertainty as large as 100% may be inferred. Then 
the fractional error in r n7f is 5r n n /r nR ~ 53%. 

We conclude that the uncertainties 5Ta/Ta and <5r nS /r n7 j on the oscillation time lower 
limits obtained in this work may be as large as 34% and 53% respectively. Of course, 
comparable uncertainties also apply to other published limits on nn oscillation times. 



8 Comparison with Previous Experiments 

Table [| compares our results with those obtained by the three most recent of previous 
searches. The numbers of candidate events and corresponding estimates for background 
rates are shown in columns four and five. Background estimates for the underground exper- 
iments (1.1, 2.5, 4.5 events for Kamiokande, Frejus, and Soudan 2 respectively) are based 
upon simulations of atmospheric neutrino interactions in the detectors. For the Frejus exper- 
iment, an alternative estimate of 2.1 background events was obtained based upon interactions 
recorded using planar spark chambers exposed to beams of neutrinos and antineutrinos at 



the CERN-PS | 40fl . Soudan 2 is the underground experiment with the largest exposure. It 
is also the only experiment to observe candidate events, although all of the underground 
experiments estimated about one background event per kiloton year. 

A degree of caution is warranted when comparing oscillation time lower limits for free 
neutrons as calculated by underground experiments versus a reactor neutron beam exper- 
iment (last column of Table |^). The r n n limits of Soudan 2, Frejus, and Kamiokande are 
based upon values for the neutron suppression factor Tr as calculated by Dover et al. [[7]. 
However if Tr values for iron and oxygen are considered which are at the high end of ranges 



obtained by Alberico et al. ||1 1|| , then limits at or below the Grenoble value are implied for 
the underground experiments. 



9 Summary And Conclusion 

The fine-grained Soudan 2 tracking calorimeter has been used in a search for neutron-to- 
antineutron oscillations occurring with bound neutrons. Based upon a fiducial exposure of 
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5.56 kiloton-years of the underground detector, a new oscillation time lower limit for nn 
oscillations in iron nuclei has been determined at 90% CL: 

T A (Fe) > 7.2 x 10 31 y. (3) 

Assuming the suppression factor for iron is Tr = 1.4 x 10 23 s _1 j7|, the corresponding 
limit at 90% CL for n to n oscillations of free neutrons is 

T nTl > 1.3 X 10 8 S. (4) 

The search reported here is background-limited. Candidate events are observed to oc- 
cur at a rate similar to that predicted for backgrounds. Since the predominant background 
arises from kinematic overlap with multiprong reactions initiated by atmospheric neutrinos, 
it seems unlikely that future underground experiments can avoid also becoming background- 
limited in larger exposures. Thus reactor neutron beam experiments, rather than under- 
ground experiments monitoring bound neutrons, may offer a more promising route for future 



improvements in sensitivity to r n ^ (see ET| and references therein). 
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